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ABSTRACT 

A graphical and an algebraic demonstration is made to show why the slope 
and zero point of the Cepheid period-luminosity (P-L) relation is rigidly coupled 
with the slope and zero point of the Cepheid instability strip in the HR diagram. 
The graphical demonstration uses an arbitrary (toy) ridge line in the instability 
strip, while the algebraic demonstration uses the pulsation equation into which 
the observed P-L relations for the Galaxy and the LMC are put to predict the 
temperature zero points and slopes of the instability strips. Agreement between 
the predicted and measured slopes in the instability strips argue that the observed 
P-L differences between the Galaxy and LMC are real. In another proof, the 
direct evidence for different P-L slopes in different galaxies is shown by comparing 
the Cepheid data in the Galaxy, the combined data in NGC3351 and NGC4321, 
in M31, LMC, SMC, IC 1613, NGC3109, and in Sextans A+B. The P-L slopes 
for the Galaxy, NGC3351, NGC4321, and M31 are nearly identical and are the 
steepest in the sample. The P-L slopes decrease monotonically with metallicity 
in the order listed, showing that the P-L relation is not the same in different 
galaxies, complicating their use in calibrating the extragalactic distance scale. 



Subject headings: stars: variables: Cepheids — P-L relations — distance scale 
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INTRODUCTION 



There is evidence that the Cepheid period-luminosity relation is not universal but dif- 
fers in slope and z ero point from galaxy-to-galaxy at a level of up t o ~ 0.3 mag as a func- 



tion of p eriod (cf . ITammann fc Reindll 120021 ; ITammann et al. 



hereafter ITSR03I : iSandage et al.ll20o3 hereafter ISTR04I 



2002 



Tammann et al 



Sandage fc Tammannl 120061 for a 



2003 



review). Drastic as this conclusion is for st udies of the extraga l actic distance scale , it ha s 
been strengthened in co n firming stud i es by iNgeow et al.l (20031). iKanbur fc Ngeowl (120041 ). 
Ngeow fc Kanbur! (booi l2005l . l2006h . INgeow et al.l (boosh . iKoen et all (l2007h . and from 
theoretical mod els a s a function of chemical co mposition by many autho r s, starting per 
haps w i th John 



Cox 



(119751). IChiosi et al 



19591 . Il980h . and includin g Ehristv et all Jl966l. Il972h. Ilben & Tugg 



and more recently iBono et al 



Fiorentino et al 



(I2002h 



Marconi et al.l (120051 ) . and undoubtedly others. 



These studies show that the position of the borders of the the L, T e instability strip in 
the HR diagram depends on chemical composition. If the strip borders vary in position and 
slope, so must the slope and zero point of the P-L relation, as worked through the pulsation 
equation in the following sections. 

Despite this evidence, the conclusion that different P-L relations apply in different galax- 
ies has recently been challenged in the literature. In these papers it is said that the slopes of 
the Cepheid P-L relations in other galaxies satisfy the slope of the P-L relation in the LMC 
and therefore that no slope differences with LMC have been demonstrated conclusively (cf. 
Gieren et al]l2005al lbl.l 20061 ; iPietrzynski et al.ll2006l ; iBenedict et al.ll2007l ; Ivan Leeuwen et al.l 
20071 are examples) . 



However, this claim sets aside the parallel evidence that the slope and zero point of the 
ridge lines of the Cepheid instability strips of the Gal axy, LMC and SMC themselves differ 
in temperature at a given period (cf. Fig. 3 of ISTR04I ). and hence, in luminosity. 



The purpose of this paper is to again remind us tha t the slope of the P-L relation is 
rigidly coupled with the slope of the instability strip via the iRitterl (118791 ) pulsation condition 
that Pyfp = constant. Hence, if the instability strip slope varies from galaxy-to-galaxy, so 
must the P-L slope. 



Differences in the instability strips of the Galaxy and SMC were first set out bylGascoigne fc Kron 

( 19651 ) . They were made secure as temperature differenc es by Lanev fc Stobie 1 jl986h . and 



have now been made definitive by the n ew CCD data by lUdalski et al.l (Il999al ibl) for LMC 
and SMC and by iBerdnikov et al.l (120001 ) for the Galaxy, as summarized for the Galaxy and 
LMC in Figure 20 of lSTR04 . 



- 3- 



In the next section we show the pulsation equation graphically and demonstrate from 
it the stated premise; a slope difference in the ridge line of the instability strip leads to a 
slope difference in the P-L relat ion. Th e grap hical sol ution here is parallel to the algebraic 

7.3: ISTR 04l S 



demonstration given elsewhere (1TSR03 



and made more explicit here in 



2. A GRAPHICAL SOLUTION BASED ON THE LINES OF CONSTANT 

PERIOD IN THE HR DIAGRAM 



In an obvious way the Ritter P^fp pulsation condition can be put into the observable 
parameters of period, luminosity, mass, and temperature by also using the Stefan-Boltzmann 
black body radiation condition that L ~ i? 2 T e 4 . The Ritter plus black body condition 
is improved by model calculations for real stars by using details of the pulsating stellar 
atmosphere structure, leading to the more precise pulsation equation of P(L, M,T e ). 



As in previous papers we again use the van Albada-Baker (119731 ) pulsation equation. 
Although it was calculated by them to apply to the lower mass RR Lyrae stars, com- 
parisons show that their predicted P-L relation is nearly identical with many other pul- 
sation equatio n s calc ulated for high er mass Cep h eids. Examples are the equations by 



Iben &: Tugglel (119751 . t heir eq. [3]). IChiosi et al.l (119921 . their eq. [5]), ISimon &: Clement 



jl993l . their eq. [2]), andlSaio fc Gautschvi Jl998h. The near identity among the equations is 
discussed in Sandage et al. ( 1999 . hereafter ISBT 99 ) . 



The pulsation equation by Ivan Albada fc Bakerl is 



logP = 0.84 log L bo i - 0.68 log Mass - 3.48 log T e + 11.502. 



(1) 



It can be made into an equation, P(L, T e ), for the lines of constant period in the HR diagram 
once a mass- luminosity relation for Cepheids is used to eliminate mass from equation ([1]). 

Observational determinations of many Cepheid masses are not available, and we must 
rely on theoretical mass values from calculated evolution tracks that pa ss throu gh the in- 
stability strip. A summary of such tracks is given in Tables 1 to 5 of ISBT 99l for tracks 
calculated from the Geneva models, in Table 11 for the Padua trac ks, and Table 12 for the 
Saio-Gautschy tracks . Detailed references for the se models a r e in ISBT 99l . The models of 
Marconi et al. ( 2005 ) for solar metallicity and by Bono et al. ( 200ol ) for lower metallicities 
were also studied. 



From all the models, normalized at log Mass = 0.84 at logL = 3.80, we have adopted 



log Mass = 0.300 logL bo i - 0.300 



(2) 
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from the tracks. This is everywhere within A log Mass = 0.03 dex of the Geneva and Padua 



tables in SBT 99 for all metallicities. 



Putting equation (j2J) into equation (TTJ gives the equation of the lines of constant period 
to be 

logL bo i = 5.472 log T e + 1.572 log P - 18.406. (3) 

This produces a family of lines in the logL, logT e HR diagram as logP is varied. 

Figure [Tj shows such a family for log P values of 0.4, 0.7, 1.0, 1.3, and 1.6. The ridge-line 
insta bility str ip for the Galaxy is shown, using its equation of logT e = —0.054 log L + 3.922 



from ISTR 04l . Figure 20. The blue and red strip borders are arbitrarily drawn parallel to the 



Galaxy ridge line using a temperature width of A logT P = 0.06 from the Galaxy ridge line. 



This is slightly wider than is observed (Fig. 20 of ISTR 04h . but is drawn to accommodate 



the dashed strip line of a toy galaxy shown with the equa tion log T P = —0.100 log L + 4.103 



similar to the instability strip of LMC (again Fig. 20 of |STR04j), but drawn here without 
the break at 10 days. The toy galaxy strip (the dashed line) has been made to intersect 
the Galaxy strip at logP = 1.3 to insure that the separate P-L relations also cross at this 
period. 

The ridge-line P-L relations are obtained in an obvious way by reading the logL (ordi- 
nate) values at the intersections of the instability strip with the constant period lines for both 
the Galaxy and the toy model. The fact that the resulting ridge-line P-L relation obtained 
for the Galaxy in this way differs from that of the toy model because of the different slopes 
of the instability strip ridge lines is obvious from this construction. 

For logL < 4.0, the instability strip has higher temperatures for the toy model than 
for the Galaxy at a given period. Hence, the intersection of the ridge-line strip with the 
constant period lines occurs at brigher lumimosities for the toy than for the Galaxy, giving 
a P-L relation for the toy model that is brighter than for the Galaxy for all periods smaller 
than logP = 1.3. The opposite is true for logP > 1.3. Hence the P-L relations will have 
different slopes, as was to be shown. 

The discussion here in words could complete the promised demonstration that the slope 
of the P-L relation is rigidly coupled with the slope of the instability strip ridge line in the 
HR diagram. However, to make the point more explicit, even to the point that the discussion 
becomes unnecessarily more elementary, bordering on pedantry, Figure [2] displays the two 
different P-L relations obtained by reading Figure [TJ in this way. The slope values for the 
Galaxy and the toy are marked in the Figure, based on the adopted instability equations 
adopted for Figure [Tj Thes e slopes are similar to the ac tual slope values measu red for the 



Galaxy and the SMC from ISTR 04j (their eq. [17]) and iTammann et al.l (120081 ). hereafter 
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TSR 08l (their eq. [5]), and set out again in Tabled] here later. 



3. THE ALGEBRAIC SOLUTION USING DATA FROM THE GALAXY 

AND THE LMC 

We can apply the pulsation equation directly to show the algebraic solution for the 
same problem using real data, both for the equations of the instability strips of the Galaxy 
and LMC and the observed P-L r elations . The demonstration made here uses the equations 



for observed P-L relations from ISTR 04l in their equation (17) for the Galaxy and their 
equations (12) and (13) for the LMC. These are put into equation ([T]), which, together with 
the adopted mass-luminosity equation (|2J), gives a predicted logT e , logL instability ridge-line 
relation. This predict ed line is then compared with the observed instability strip equations 



shown in Figure 20 of ISTR04 . 



We have used an explicit bolometric correction to change the \ogL v values obtained 
from the observations into logLboi required in equations and (jSJ), and back to logLy to 
compare the predictions from the pulsation e quation w ith the observations. The bolometric 



corrections are interpolated from Table 6 of iSBT 99l for the appropriate metallicities and 
surface gravities of the Cepheids. The turbulent velocity was assumed to be 1.7km s - . 
The surface gravities vary with radius, mass, and l uminosit y and therefore with period as a 



surrogate as logg = — 1.09 log P + 2.64 (eq. [49] of [STROJ). The metallicities are assumed 
to be [A/H]= 0.00 for the Galaxy and —0.5 for LMC. The mass is from equation (121) . The 
obvious arithmetic is not shown. 

The resulting predictions of the instability strip ridge-line equations are these: 

log T e (predicted) = -0.040 log L v + 3.854 (4) 

for the Galaxy at all periods, and, 

log T e (predicted) = -0.056 log L v + 3.941 (5) 

for P < 10 days for the LMC, and 

log T e (predicted) = -0.081 logLy + 4.020 (6) 

for P > 10 days, also for the LMC. 

Note that the break in the T e —L instability strip relation at P = 10 days in equations flSJ) 
a nd (El) is mirrored in the break in the P-L LMC relations given in equations (12) and (13) 
of ISTR04 , and shown as Figure 4 there. 
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For comparison with the predictions in equations (jl])-© here, the observed ridge lines 
of the strip s in the Galaxy and the LMC, taken from the insert equations shown in Figure 20 
of lSTRpj are 

log T e (observed) = -0.054 log L v + 3.922 (7) 
for the Galaxy at all periods, and 



for P < 10 days, and 



log T e (observed) = -0.050 log L v + 3.936 



log T e (observed) = -0.078 log L v + 4.029 



(8) 



(9) 



for P > 10 days for the LMC. 

The near agreement of the predicted slopes of the instability strips in equations (JI])-© 
with the observed slopes in equations dTJ)-© is the demonstration we are seeking. 

The agreement is good, but there is a disagreement in the temperature zero points 
between equations PJ-([S5 and equations ([T])-® by AlogT e = 0.018 dex. The predicted 
temperatures are cooler than those observed. However, the difference is remarkably small, 
given the approxi mations we have made in the bolometric corrections, in the adopted tem- 
perature scale of ISBT 99i . their Table 6, and in the adopted van Albada-Baker theoretical 
zero point in equation (jTJ). 

The temperature offset could be made zero if the zero point of the mass in equation ([21) 
would be made smaller by 0.09 dex, but then the evolution mass would differ from the 
pulsation mass by this amount. This is the expression of the previous well known mass 
"problem" which is solved here by the temperature shift. 

In this regard, it is useful to remark t hat many of the temp e rature scales in the current 



l iterat ure, for example as summ arized b y ISekiguchi fc Fukugital (120001 ) or by ICacciari et al. 



(120051 ). and including the one in ISBT 99l that we have used here, differ among themselves by 
as much as 0.025 dex in logT e at fixed B — V. Th is, then, is the temperature uncertainty in 
the temperature zero point in Figure 20 of |STR04|. Our shifting of the predicted temperature 
relative to the observed temperatures in Figure [3] by 0.018 dex is not excessive. 

The observed (solid lines) and the predicted (dashed lines shifted by 0.018 dex in logT e ) 
instability strips for the Galaxy and the LMC are shown in Figure [3J The agreement is 
satisfactory, showing again that differences in the instability strip loci causes differences in 
the slopes of the P-L relations. Hence the claims in the current literature, cited in the 
Introduction, that a universal slope exists for the Cepheid P-L relation are inconsistent with 
Figures [TH3] which show different positions of the instability strip in different galaxies. 
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SUMMARY OF OBSERVED P-L SLOPE DIFFERENCES IN SELECTED 

GALAXIES 



The arguments given in the previous sections rely on knowledge of the temperatures 
of the instability strips. These can only be measured using reddening corrected colors, and 
these are reliable only if the reddening of the individual Cepheids can be determined by some 
method other than by using a fiducial period-color (P-C) relation. The reason is that if the 
temperatures of the instability strips differ from galaxy-to-galaxy, presumably because of 
chemical composition differences, the P-C relations will also differ. There will be no correct 
fiducial P-C template from which to determine the reddening if the chemical compositions 
vary greatly. The intrinsic P-C relations will differ from galaxy-to-galaxy depending on the 
chemical composition, and the reddenings are therefore indeterminate. 

Presently, it is only the Galaxy, LMC, and SMC that can be subjected to the analysis 
given here because it is only for these galaxies that the reddening of their Cepheids have been 
determined by methods other than by comparing with some adopted fiducial P-C relation. 

However, for some galaxies with enough Cepheids, and where the differential reddening 
between the Cepheids is small enough to be ignored, comparison of the P-L slopes can be 
made directly from the data. The result for the Galaxy, NGC3351, NGC4321, LMC, SMC, 
IC 1613, NGC 3109, and Sextans A and B is shown in FigurelU The adopted data for the P-L 
relations are in Table [H The equations for the apparent magnitude and absolute magnitude 
P-L relations are V° = alogP + b, and M v = a logP + c. Column (2) is the log of the 



oxygen-to-hydrogen ratio from Table 4 of |TSR 08l . Column (3) shows the observed value of 
a, which is the slope of the apparent magnitude P-L relation taken fr om the same sources 
(but changed slightly in a few entries here) as were used for Table 4 of ITSR08I . Column (4) 
lists the apparent magnitud e P-L int ercept, b, as observed. Column (5) lists the (m — M)° 
distance modulus adopted in |TSR 08l . The absolute magnitude P-L relation is in column (6), 
which is column (4) minus column (5). The literature source is in column (7). The resulting 
P-L relations, calculated from the a and c values in Table [TJ are shown in Figure HI The 
slopes for the NGC 3351/ NGC 4321 combination and the Galaxy are the steepest of those 
shown, and are similar. That of the LMC is next steepest. 



The slope of the Galaxy P-L relation in |TSR 03l and ISTR 04 is based on averaging the 
results using the moving atmosphere method (the Baade-Becker-Wesselink procedure) and 
the independent main sequence fitting method. Ne vertheless, the result i ng slope of the P-L 



slope has been questioned as being too steep (cf. iGieren et al.l l2005bl ; Ivan Leeuwen et al. 
20071 ). Howeve r, the slopes of the N GC 3351 and NGC 4321 combined P-L relation, and 



th at of M31 bvlVilardell et al.l (120071 ) are eq ually ste ep as for the Galaxy. The M31 slope 
by IVilardell et al.l has been redetermined by ITSR08I . The original slope by Villardel and 
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collaborators was based on E(B — V) values using the LMC P-C relation rather than the 
more correct higher metallicity P-C relation for the Galaxy. The resulting E(B — V) values 
turns out to depend on period as a further complicati on. But even discounting the M31 
case, the steep slope for NGC 33 51 and NGC 4321 from lTSR08l (their Fig. 2), supports the 
Galaxy slope that we derived in ISTR 04 and its difference from the P-L slope in LMC. 



The s trongest evidence for t he difference as function of metallicity is the data for 
NGC 3109 Jpietrzvnski et alil2006h which has the well determined P-L slope of dM v /d log P = 
—2.13. This differs significantly from the slopes of either the Galaxy or NGC 3351/ NGC 4321 
at —3.10, or —2.92 for M31, and even for the LMC at —2.70. The six longest period Cepheids 
in NGC 3109 with logP > 1.3 are too faint by ~0.2 mag (Fig. 4 of IPietrzvnski et alJliooS ) 
compared with either the Galaxy or the LMC P-L relations. 



Figure HI similar in principle to Figure 5 in |TSR 08l . together with Figure [3] here, is our 
chief case for non-unique P-L relations between galaxies of different chemical compositions. 
The complications that this portends for determining the scale of extragalactic distances 
from Cepheids to withi n ~ 15%, unless special corrections for the diff erence are applied, is 
discussed elsewhere (cf. ISaha et al.ll2006l ; ISandage et al.ll2006l ; ITSR081 ). 



It is a pleasure to thank Bernd Reindl for his skill in the preparation of the diagrams 
and the text for publication. We also thank John Grula, Carnegie editorial chief, for his 
liaison with the press. 
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Table I. Observed (P-L)y Relations for Ten Galaxies with Different 

Chemical Compositions 



Name 


[O/H] 


a 


b 


(m - Mf 


c 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) ' 


(6) 


(7) 


Galaxy 


8.60 


-3.087 






-0.91 


1 


NGC 3351/4321 


(8.80) 


-3.108 


mean 


mean 


-0.90 


2 


M31 


8.66 


-2.92 




24.43 




3 


LMC 


8.34 


-2.702 


17.05 


18.54 


-1.49 


4 


SMC 


7.98 


-2.588 


17.53 


18.93 


-1.40 


5 


IC1613 


7.86 


-2.698 


23.08 


24.35 


-1.27 


6 


NGC 3109 


8.06 


-2.130 


23.73 


25.45 


-1.72 


7 


Sextans A/B 


7.52 


-1.628 


23.10 


25.80 


-2.40 


8 



References; — ( 1 ) [STRoJ, eg. [171; f2llTSR08L Fie. 2; ( 31 IVilardell et al 



2007: (4) STR04 


eq. 


Pietrzvnski et al. 


2006 



8]; f5l ITSR 08l. ea. [51: ( 61 



Antonello et al 



20061 : (7) 



4; (8) lPiotto et alJll994 . 
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Fig. 1. — Schematic HR diagram in the vicinity of the Cepheid instability strip. The central 
line is the observed ridge line for the Galaxy taken from Figure 20 of STR 04 whose equation 
is logT e = —0.054 log Ly + 3.922. The dashed line is for a toy galaxy whose ridge-line 
equation is logT e = — 0.100 log Ly + 4.103. The borders of the instability strip are put 
parallel to the Galaxy ridge line. Lines of constant period, calculated from equation ([3]), are 
marked with their logP values (in days). 
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Fig. 2. — The two P-L relations for the two ridge lines in Figure HJ determined from the 
intersections of the ridge lines of the Galaxy and the toy galaxy with the lines of constant 
period in Figure [TJ The absolute magnitudes along the ordinate are transferred from Figured] 
by My = — 2.51ogL bo i + 4.75 where the bolometric correction in V is adopted to be zero. 
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Fig. 3. — The algebraic demonstration of the rigid coupling between the slopes of the insta- 
bility strip and the slope of the P-L relation required by the pulsation equation. Predicted 
(dashed lines) slopes and zero points for these instability strip ridge lines i n the G alaxy 
and the LMC are compared with the observed (solid) lines from Figure 20 of ISTR 04 The 
predictions are made by inserting the equations of the observed P-L relations for the Galaxy 
and the LMC into the pulsation equation flTJ. The predicted zero points are moved by 0.018 
in logT e , hotter. 
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Fig. 4. — The observed ridge lines of the P-L relations for eight galaxies listed in Table [TJ 
The P-L relation for the Galaxy (not shown) is nearly identical with the combined NGC 3351 
and NGC 4321 line, and has the steepest slope. The agreement between the Galaxy and the 
combined NGC 3351 and NGC 4321 slopes argues for the correctness of the steep slope for 
the Galaxy P-L relation. 



